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Abstract 

Equations of motion for a classical 3d discrete model, whose auxialiary system is a linear system, are 
investigated. The Lagrangian form of the equations of motion is derived. The Lagrangian variables are 
a triplet of "tau functions" . The equations of motion for the Triplet of Tau functions are Three Trilinear 
equations. Simple solitons for the trilinear equations are given. Both the dispersion relation and the phase 
shift reflect the triplet structure of equations. 

1 Introduction 

Recently there was proved the exact integrability of one discrete 3d model, classical 
as well as quantum [Q, ^. This model was formulated originally in terms of the 
pure time evolution as the map of the dynamical variables from time t to time t + 1. 
Namely, for t fixed the system of the dynamical variables are the system of 

[^a,a,b , ^a,a,b] , « = 1, 2, 3, a,b £ Zm , (1) 

where M is the finite spatial size of the system, and the Poisson brackets in the 
classical case are 

{Ua,a,b , ^a',a',b'} = ^a,a' Sa,a' ^b,b' ^a,a,b^a,a,b ■ (2) 

The evolution was formulated as the explicit form of functions /, g, 

[u,w] = [u(t),w(t)] ^ [u(t+l),w(t + l)] = [/(u,w),5(u,w)] , (3) 

such that the Poisson brackets are conserved by this map. The existence of the 
complete set of involutive integrals of motion was proved in |l|, • 

The map t i— s- t + 1 may be considered as a sort of the Hamiltonian equations of 
motion for this classical discrete model. In this paper we will derive their Lagrangian 
counterpart. Following one may guess a priori that the Lagrangian variable is 
Hirota's tau function and it might obey Hirota's discrete bilinear equations 
But we will show that this is not true in general, although equations we will derive 
resemble Hirota's ones. 

In this note we recall the reader the Hamiltonian form of the equations of motion 
first, then introduce the analogous of the tau function and derive the Lagrangian 
form of the equations of motion, and finally give the multi-soliton solution of them. 



2 Hamiltonian equations of motion 



In this paper we investigate the equations of motion, so we need no the evolution 
operator and may choose more appropriate coordinate system. Actuahy, the co- 
ordinates we will use here are the light cone frame with respect to previous t, a, b. 
Let 61,62,63 be three translations making 



61 : {t,a,b) 1-^ (t + l,a, 6) , 

62 : {t,a,b) {t + l,a + l,b) , 

63 : (t,o, 6) 1-^ (t + l,a, 6 + 1) . 



(4) 



61 ) 62 ) 63 are the elementary orthonormal shifts of the three dimensional cubic lat- 
tice. So for any site of the cubic lattice p given, the following eight sites form the 
elementary cube of the lattice: 



^ P , p -I- 62 63 , p 61 63 , p -I- 61 62 , 

y p 61 62 + 63 , P + 61 , P + 62 , P + 63 . 



(5) 



Equations of motion for Ua,p, Wq,,p may be extracted from the form of f,g in 
the definition of the evolution in |]l|, §], and in our coordinates look like 
(i) 

Wl,pW2,p + U3,pW2,p + K3U3,pW3_p 
Wl,p+6i = 



Ul,p+6i 



W3,p 
fi2 Ul,p U2,p W2,p 



(6) 



(ii) 



W2,p+62 
U2,p+62 



(iii) 



W3,p+63 



^J3,p+63 



KlUi_pW2,p + K3U2,pW3_p + Ki K3 Ui^p W3^p 

_ Wl,pW2,pW3^p 

Wl,pW2,p + U3pW2,p + K3U3^pW3p 

^ Ul.p U2,p US.p 

U2,pU3,p + U2,pWi,p + KlUi,pWip ' 

K2 U2,pW2,p W3p 

KlUi,pW2,p + K3U2,pW3,p + Kl ^3 Ul,p W3,p 

U2,pU3,p + U2,pWi,p + KiUi^pWi^p 

ui,p 



(7) 



(8) 



Describe in a couple of words the idea of the derivation of equations These 
relations may be obtained as the zero curvature condition of the following system. 
Let ipp be an auxiliary variable assigned to the sites of the cubic lattice. Consider 
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three orthogonal plaquettes of the cube endowed by the following relations 
U = Ji,p = wi^p - (/9p+e2+e3 ui,p + ^^p+ea + f/^p+ea i^i ui,p wi^p , 

U = /2,p = y'p + (/'p+ei+es ^^2 U2,p W2,p + y'p+ei W2,p - y'p+ea U2,p , 

def 



= /3,F 



V'p U3,p + 99p+ei+e2 W3,p + 93p+e2 + V'p+ei ^3 U3p W3p . 



(9) 

Easy to see that the number of linear equations is three times greater then the 
number of (p^. This means that the coefficients of the system of linear equations 
must obey extra conditions. As an example consider a cube. Six plaquettes give 
six linear relations for eight corner </?p: 



/o,p — /a,p+eQ — 0, a — 1,2,3 



(10) 



Demand that between these six relations there are only four linearly independent, 
so that the cube relations fix only four of eight corner linear variables (/?p (e.g. we 
just may express (/^p+ei+ea, Vp+ez+es, V^p+ei+es and (/Jp+ei+ea+ea via independent 
93p, 93p+eii ^p+e2 and v^p+es and nothing more). This demand gives equations of 
motion (|6[|1) immediately. 



3 Lagrangian equations of motion 

Three simple relations between the equations of motion are to be mentioned: 

Wl,pW2,p = Wi^p+ei W2,p+e2 ; 

U2,p U3,p = U2,p+e2 U3,p+e3 ) (11) 

. ui,p/w3^p = ui_p+ei/w3,p+e3 . 

Following 1^ , a parametrization of u , w in terms of tau functions must turn relations 
(|Tl| ) into tautologies. Thus without lost of generality (we deal with the infinite 
system, avoiding hence all the boundary problems) 



w 



i>P 



w 



2,P - T, 



w 





ui,p = 


^,P 


'^3,p ' 


'^.p+es 


^,P 


U2,p = 


^i,P 


^,p+ei 


'^i.p+es 


'?2,p 


^3,p = 


^i,p+e2 


'^.p+ei 


^i,P 



(12) 



3,P - ^ 

As usual, 7^,p may contain a quadratic and linear pre-exponent, 

r n = e5(P'«"'P) + (qc.P) . r n 



(13) 
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where q^, and diagonals of Qa are arbitrary, but all non-diagonal elements of Qa 
must coincide (in the natural basis e^). 

Substituting now the parametrization (|l^) to (|6[]^ and taking ( p^ ) into account, 
we obtain Three Trilinear relations for Three r functions: 



ri ri,p+e2+e3 ^2,^ r3,p 



^2 Ti,pr2,p+ei+e3 T3,p 



Ti,pr2,pr3,p+ei+e2 



= Ti,p T2,p+e3 "^a.p+ea 

+ S2 ''"i,p+e2 'r2,p+e3 '^3,p 

+ S3 ^ n,p+e3 'r2,p T3,p+e2 

= ''"i,p+e3 T2,pT3,p+ei 

+ S3 Ti^p r2,p+e3 T3^p+ei 

+ s^^ n,p+e3 'r2,p+ei ■^s.p , 

= n,p+e2 T"2,p+ei T3,p 

+ Sx ri^p+e2 ''"2,p''"3,p+ei 



(14) 



(15) 



(16) 



+ S2 Tipr2,p+ei T-3,p+e2 • 
Here we gather all k-s and all the pre-exponents into arbitrary Sa,rQ,, a = 1,2,3. 
One may change them in any appropriate way. 

Because of the possibility to introduce the Poisson brackets (^) in the system of 
Uo,p, Wa^p, the original equations of motion (§-§) may be regarded as a kind of the 
Hamiltonian equations of motion. Contrary to that, the tautological parametriza- 
tion of ( [11] ) resembles the introduction of a velocity instead of a momentum, so 



that the number of variables decreases twice, so the trilinear relations (14-1^ are 
nothing but the Lagrangian form of the equations of motion. 



(14-16D may be reduced to Hirota's bilinear discrete equation in the limit 

Kl << K2 = «^3 << 1 • (17) 



4 Solitons 

Now describe the simple solitons of the trilinear relations (|T^^). 
Let (a,/?, 7) be any cyclic permutation of (1,2,3). Fix via 

= I + sp + ^ . (18) 

Then let the exponent 

^ = e^^'^'P) = e^^i^i +^^2P2 + ^^3P3 (19) 
for Sa given, is defined by 

A/3 (Aq — A.y) + Sa (Aq, — Xfs) 
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This parametrization of the dispersion relation we will exhibit as W = W{X). For 
the dispersion curve W = W{X) given, let for the shortness 



The phase shift D is defined by 

^(A,A') 

where 



d{X,X') ■ d{X-\X'-^) 



d{X-\X') ■ (i(A,A'-i) 
/ 1 1 1 \ 



d(A,A') 



det 



Ai Ao Aq 



V A'l A2 A3 / 

and A~^ stands for {A^^, A^^, A3 ^}. For the shortness let 
One soliton solution is rj^p = 1 — Wa- Two solitons are 

.(2) 



Three solitons are 



J3) 



+ Z)(1'2) T^^l) + ^(1,3) ^(1) ^(3) ^ ^(2,3) ^(2) ^(3) 

- D(1'2) Z)(1.3) £,(2,3) wW wi^) . 



(21) 



(22) 



(23) 



(24) 



(25) 



(26) 



And so on, common expression is usual. Remarkable is that the dispersion rela- 
tion for the exponents e*^" may be parametrized by their origin A^, (pO|) , and 
parameters Sa are defined by an arbitrary linear pre-exponent. 



5 Discussion 

System (p"4|-[T6D contains a set of bilinear relations as its necessary conditions. Look- 
ing for a solution of trilinear relations in a form of holomorphic functions, we get 
the following set of necessary relations 

Aa.p+e^, T-Q.p = Ti3,p+ea^j,p + •5a T/jp r^,p+e^ , (27) 

where as before (a, (3, 7) are any cyclic permutation of (1, 2, 3), and auxiliary func- 
tions Aa^p are supposed to be holomorphic. Thus (ll^p^) have a host of solutions 
similar to that of Hirota's equation. Namely, identifying (27) with a set of Fay's 



or Pliicker's relations, we may obtain elliptic or determinant solution of (jT^|T^). 
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Probably, the exception is a Bethe - ansatz - type solution of (p7|), because in this 



case (14 1^ do not follow from ( P7|) tautologically. Because of all these do not 
contain anything so surprising as the solitons, we do not give any explicit formula 
here. Mention only two aspects. 

First one is the symmetry group of ([T^|l^). For Hirota's equation all the sym- 
metries with respect to permutations and reflections of the space directions are 
obvious. In the case of (H-IH) the Cube group is less trivial because of it acts on 
the parameters Sq , . As an example give the realisation of the complete reflection 
P. For a solution r^p of the trilinear equations with some given let 

- Sg _ s/3 rpr^ 
Ta,p = Ta,p_e^ , Sa = , = . (28) 

Then the set of "over lined" objects obey the same set of the trilinear relations 
(|l|-|l6l) but with reflected Sq, I — > — Sq: for all OL. 

Mention also the possibility to construct the complete solution of original finite 
evolution model in terms of the elliptic functions on a finite genus curve. In such 
case one may consider a non-homogeneous and non-equidistant version of the tri- 
linear relations with a boundary conditions taken into account. The number of the 
parameters (moduli, initial divisor and homogeneousless data) will coincide with 
the number of variables of the initial state. This will be the subject of a separate 
paper. 
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